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Abstract 

We show how the pre-exponential factor of the Feynman propagator for the harmonic 
oscillator can be computed by the generalized (^-function method. Besides, we establish a 
direct equivalence between this method and Schwinger's propertime method. 
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In a recent paper that appeared in this journal [1] the harmonic oscillator propagator was 
evaluated by a variety of ways, all of them based on path integrals. In fact, some of them did 
not involve any explicit computation of the Feynman path integral, but their common starting 
point was actually an expression for the harmonic oscillator propagator which was explicitly 
derived by path integral means, namely (we are using as much as possible the notation of 
reference [1]): 
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and the determinants must be computed with Dirichlet boundary conditions. In the previous 
equation S[z d ] means the classical action, that is, the functional action evaluated at the classical 
solution satisfying the Feynman conditions z(ti) = Zi and z(tf) = Zf and the factor before the 
exponential is the so called pre-exponential factor, which we shall denote by F(tf — ti). In 
reference [1] three distinct methods were presented for the computation of F(tf — U): (i) it 
was computed directly by the products of the corresponding eigenvalues of O and (some 
care must be taken here since both products are infinite, but their ratio is finite); (ii) it was 
computed with the aid of Schwinger's propertime method [2] (an introductory presentation of 
this method with simple applications can be found in reference [3]); (in) it was computed by 
the Green function approach (a variety of simple examples worked out with this approach can 
be found in references [4, 5]). 

In this note we just add to the previous list one more method for computing the pre- 
exponential factor of the harmonic oscillator propagator, namely, the generalized ( function 



method, so that this note can be considered as a small complement of Holstein's paper [1]. 
In fact, every time we make a semiclassical approximation, no matter it is in the context of 
quantum mechanics or quantum field theory, we will get involved with the computation of 
a determinant of a differential operator with some boundary conditions. If we try naively 
to compute these determinants as the products of the corresponding eigenvalues we will get 
ill defined expressions. Hence, it is imperative to give a finite prescription for computing 
determinants for these cases. The generalized (^-function method is precisely one possible way 
of doing that. It was introduced in physics in the middle seventies [6] and it is in fact a very 
powerful regularization prescription which has applications in several branches of physics (a 
detailed discussion can be found in reference [7]). This method, as we will see, is based on an 
analytical extension in the complex plane. We think that the harmonic oscillator propagator is 
the perfect scenario for introducing such an important method, because undergraduate students 
are all familiarized with the quantum harmonic oscillators and besides, it is the first non-trivial 
example after the free particle. In what follows, we shall first introduce briefly the ^-function 
method, then we shall apply it to compute F(tf — ti) for the harmonic oscillator propagator and 
finally, we shall establish a direct equivalence between this method and Schwinger's propertime 
method. 

Consider an operator A and let us assume, without loss of generality, that it has a discrete 
set of non-degenerate eigenvalues {A„}. When there is only a finite number of eigenvalues det A 
is just given by the product of these eigenvalues and we can write: 
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where we defined the generalized zeta function associated to the operator A as 



((s;A) = TrA 



—s 
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However, when there is an infinite number of eigenvalues (and these are the cases of interest 
in physics), as it occurs when A is a differential operator, the product of the eigenvalues will 
be an ill defined quantity and will not serve anymore as a good prescription for det A. In other 
words, expression (3) with ((s; A) given by (4), as it stands, is meaningless because it is not 
valid anymore to use that: 



Hence, we need for these cases to define a finite prescription for detA. The generalized zeta 
function prescription consists basically of the following three steps: (i) we first compute the 
eigenvalues of A subjected to the appropriate boundary conditions and then write down the 
corresponding (^-function £(s; A) = Tr A~ s = J2n ^n S ! 00 Since the last sum does not converge 
at s = 0, we make an analytical extension of this function to the whole complex plane of s (or 
at least to a domain that contains the origin); (iii) after the analytical extension is made, we 
just write det A = exp {— ('(s = 0; A)}. 




(5) 
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In order to apply the ^-function method described above in the computation of F(tf — U) 
we first need to find the eigenvalues of O with Dirichlet boundary conditions. For convenience, 
we shall make the rotation in the complex plane t = e~ J7r / 2 T = —iT. Let us also define the 
corresponding finite interval in T by tf — U = —i(Tf — Tj) = —i r. We then have that 

This analytical extension guarantees that all the eigenvalues (now of the operator Or) are 
positive. Of course, after the calculations are finished we must undo this transformation, that 
is, we must substitute r = i(tf —ti). Solving the eigenvalue equation Ot f n (T) = X n f n (T) with 
Dirichlet boundary conditions / n (0) = = / n (r) we get: 

2 2 

/ n (T) = {sin (^r) ; n = l,2,...} ; A n = u 2 + ^ ; n = 1, 2, ... (7) 
Consequently, the associated generalized ^-function is given by 

where we defined z/ = ut/ti. Since the above series converges only for Ties > 1/2, we need to 
make an analytical extension in the complex plane of s to include the origin. However, this can 
be done with no effort at all, for this series is precisely the so called non-homogeneous Epstein 
function, which we shall denote simply by Ei(s; u 2 ) and whose analytical extension to the whole 
complex plane is well known and is given by [8, 7] (see the appendix for a brief deduction): 
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where K^z) is the modified Bessel function of second kind. This is a meromorphic function in 
the whole complex plane with simple poles at s = 1/2, —1/2, —3/2, so that we can take 



its derivative at s = without any problem. Substituting the sum appearing on the r.h.s. of 
(8) by the analytical extension given by (9), we may cast £(s; O t ) into the form: 



where 
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is analytic at s = 0. Taking, then, the derivative with respect to s at s = and using that 
r(s) ~ 1/s for s — > 0, we get 

f r'fsl F'(s) 1 

C'( S = 0;O T ) = -log (t/tti/) + lim \--±-Lf( 8 ) + 



s^o { T 2 (s) v ' T(s) J 
- log (t/tti/) + F(0) . (11) 



From the above expression for F(s) we readly compute F(0), so that 

C(s = 0; O t ) = - log (t/tti/) + [ v^gzl/g) + 2^ £ J^K_ 1/2 {2 
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Using that T(— 1/2) = —2^ and fT_i/ 2 (-2) = \Jn/2ze z , we obtain 
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It is not a difficult task to show that the above sum is given by (take its derivative with respect 
to is, sum the resultant geometric series and then integrate in v, in order to eliminate the 
arbitrary integration constant, just use the fact that this sum must vanish for v — > oo ): 
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From equations (13) and (14) we then have 
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where we used that v = lut/h. For the operator we immediatly get (it suffices to make 
uj — > in the previous formula): 
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Collecting all the previous results and rotating back to the real time (r = i(tf — U)), we finally 
obtain 



F(t f - U) 



exp [-C'(0,O) 



x 



m 



(17) 



where we used that sinh(i^) = —i sin 9, in perfect agreement with [1]. 

Before we finish this note, we think it is interesting to establish a general equivalence between 
the ^-function method and the Schwinger's propertime method. From the ^-function method 



just presented, we can write 



logdet O = -('(s = 0; O) . 



(18) 



On the other hand, with the aid of the Mellin transform [9] applied here to an operator O with 
positive eigenvalues we can write 



TrO- s 



= Tr 



T(s) 



Jo 



(19) 



However, last expression is not analytic at s = (though the presence of the exponential 
guarantees a well behaviour for large s, the limit s — * is a divergent one), so that as it stands 
it is not valid to take the s derivative at s = 0. In order to circumvent this problem, we make 



the modification (regularization) 
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where a is big enough to ensure that the previous expression is well behaved at s = 0. Hence, 
taking first the s derivative at s = and then taking the limit «-^0we obtain: 



Last expression corresponds precisely to Schwinger's formula written in a regularized way. Here 
we regularized by introducing positive powers of r, but other regularization schemes can also 
be used, as for example the one used by Schwinger [10] in the computation of the Casimir effect 
[11] (for a simple introduction to this effect with some historical remarks see reference [12]). It 
is common to write the above expression formally with a — 0, but in fact, before taking this 
limit one should get rid off all spurious terms (those with no physical meaning). 

In this note we have presented the generalized C-function method for computing determi- 
nants in a very introductory level. A detailed discussion with a great variety of examples can 
be found in reference [7]. One of the greatest advantadges of this method is that for almost all 
differential operators and boundary conditions that are relevant in physics, the corresponding 
generalized (^-function (after the analytical extension is made) is a meromorphic function in the 
whole complex plane which is analytic at the origin. Furthermore, this method can also be ap- 
plied successfully in many other branches of physics, as for example, statistical mechanics and 
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quantum field theory among others. Of course there are many easier ways of obtaining F(tf—ti) 
for the harmonic oscillator, but our purpose here was to introduce a new method, which is a 
powerful one and widely used nowadays. In this sense, we think that the harmonic oscillator 
provided a perfect scenario for the understanding of the three basic steps of the method, since 
every undergraduate student is somehow familiarized with the harmonic oscillator. 



Appendix 

In this appendix we shall obtain the analytical extension of the Epstein function E" 2 (s;l), 
given in the text by equation (9). With this goal, we first write down an equation involving 
the Gamma function, which follows directly from its definition, namely [9]: 

roo 

T(z)A- z = / dr r 2 - 1 e~ Ar , Re(z) > (22) 
Jo 

Using the previous equation with A = n 2 + u 2 , the Epstein function can be written in the form: 

OO 1 

Ef(s;l) = £ 



1 X (n 2 + v 2 ) 
1 

" r( 

On the other hand, from the so called Poisson summation rule [9], we can write that: 



1 r oo < ^~ > 

— / dr r 3 - 1 e- vlr V e~ n2T . (23) 
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E^ = -J + ^ + J;Ee-" V(lW (24) 

71=1 ^ ^ V T V T 71=1 



Substituting the last into (23), we get: 



Ef (s; 1) = — L [°° dr T s - x e- v2r + [°° dr r s ' 3 / 2 e~^ T 

1 V ' ; 2F(s) Jo 2T(s) Jo 



+ #E f 'dr^e-^ 1 ^ (25) 

^ \ S ) 71=1 J® 
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Using (22), the first and second integrals of the right hand side of the above equation can 
be written directly in term of Euler Gamma functions. For the last term, we use the integral 
representation for the modified Bessel function of second kind: 

J™ dxx^x-M*-^ = 2(^y /2 K a (2^hf) , Re(3,Re^>0. (26) 
Therefore, we finally obtain for equation (25) that: 

K (s: 1) = + \ T l )v Ji + W)U~> K -^ m) ' (27) 

which is precisely equation (9). Some comments are in order here: (i) to say that the previous 
equation corresponds to the analytical extension of E p2 (s; 1) to a meromorphic function in the 
whole complex plane means that this expression is an analytical function in the whole complex 
plane except by an enumerable number of poles (which can be infinite) and coincides with the 
original sum in the region where the sum was defined; (ii) it is worth emphasizing that the 
above expression is analytic at the origin; in fact, the structure of simple poles of this function 
is dictated by the poles of the Euler gamma function. It is easy to see that the poles are located 
at s = 1/2, -1/2, -3/2, -5/2. etc.. 
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